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Anomalous dispersion cannot occur in a transparent passive medium where electromagnetic ra- 
diation is being absorbed at all frequencies, as pointed out by Landau and Lifshitz. Here we show, 
both theoretically and experimentally, that transparent linear anomalous dispersion can occur when 
a gain doublet is present. Therefore, a superluminal light pulse propagation can be observed even 
at a negative group velocity through a transparent medium with almost no pulse distortion. Con- 
sequently, a negative transit time is experimentally observed resulting in the peak of the incident 
light pulse to exit the medium even before entering it. This counterintuitive effect is a direct result 
of the rephasing process owing to the wave nature of light and is not at odds with either causality 
or Einstein's theory of special relativity. 



PACS number(s): 03.65.Sq, 42.25.-p, 42.25.Hz, 42.50.-p 



I. INTRODUCTION 

In a seminal paper, Lord Rayleigh remarked that a 
pulse of light travels at the "group velocity" instead of 
the phase velocity inside a medium [Q. In subsequent 
papers, Lord Rayleigh further developed the theory on 
transparency and opacity, and the theory of anomalous 
dispersion Anomalous dispersion was first stud- 

ied for mechanical oscillators [3] and was later applied 
by Sommerfeld and Brillouin [Q to light propagating in 
absorptive opaque materials. They showed theoretically 
that inside an absorption line, the dispersion is anoma- 
lous, resulting in a group velocity faster than c, the vac- 
uum speed of light. Such an anomalous velocity appears 
owing to the wave nature of light |||| ■ ■ 

This however raised difficulties with the theory of rela- 
tivity which states that no velocity can be higher than c, 
the velocity of light in a vacuum. Brillouin pointed out 
Q "Group velocity, as originally defined, became larger 
than c or even negative within an absorption band. Such 
a contradiction had to be resolved and was extensively 
discussed in many meetings about 1910." 

In order to reconcile the superluminal group velocity in 
an anomalous dispersion medium with the implied limi- 
tation from relativity, Sommerfeld and Brillouin pointed 
out that causality only requires the speed of a signal car- 
ried by light be limited by c, rather than the light pulse 
itself which travels at the group velocity Q| . They further 
pointed out that the correct definition of the speed of a 
light signal should be defined as the "frontal velocity," 
instead of the somewhat misleadingly named "signal ve- 
locity" defined as the speed of the half-point of the front 
edge of a light pulse. A "frontal velocity" marks the ve- 
locity of an infinitely sharp step-function-like change in 



the light intensity, albeit even in principle such an in- 
finitely sharp change in time requires infinite bandwidth 
and thus becomes impractical. It was shown that since 
the infinitely sharp "front" contains infinitely large band- 
width, a good portion of the power would exceed any 
practical resonant and non-resonant (plasma) frequen- 
cies to have an effective refractive index of unity. Thus 
this portion of the power will propagate exactly at c to 
exit the medium at the earliest time and form "precur- 
sors" Q| . It was an ingenious argument that resolved the 
apparent conflict of special relativity and superluminal 
group velocities. 

More recently, it has been considered |@hTi],[13] [l|,|2l] Q 
to use various schemes including the anomalous disper- 
sion near an absorption lin e [|s| - [ll"l| , nonlinear [fl3f and 
linear gain lines li~i|- |l9| , |2^ -|24|| , active plasma medium 
|p7| , pT| , or in a tunneling barrier P^j^Jto observe these 
effects. Sometimes it is applied 2f]]2(| to observe 



superluminal propagation of light pulses. Inside an ab- 
sorption line, it was shown that a light pulse propagates 
at group velocities faster than c and can become negative 
with dramatic values of -c/23, 000 § [[J. However, in 
all experiments light pulses experienced either very large 
absorption p,^0| or severe reshaping |25|] that sometimes 
resulted in controversies over the interpretations. 

In a series of papers |]7|, ^8| , ^9|j2^ -|2rj] , Chiao and cowork- 
ers showed theoretically that anomalous dispersion can 
occur inside a transparent material, particularly on the 
DC side of a single gain line. It was predicted that by us- 
ing a gain doublet |22j] , it is possible to obtain a transpar- 
ent anomalous dispersion region where the group velocity 
of a light pulse exceeds c with almost no pulse distortion. 

Here we use gain-assisted linear anomalous disper- 
sion to demonstrate superluminal light pulse propaga- 
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tion with a negative group velocity through a transpar- 
ent atomic medium [ p7| . We place two Raman gain peaks 
closely to obtain an essentially lossless anomalous disper- 
sion region that results in a superluminal group velocity. 
The group velocity of a pulse in this region exceeds c and 
can even become negative [^2|,^3| , while the shape of the 
pulse is preserved. We measured a negative group ve- 
locity index of n g — — 315(±5). Experimentally, a light 
pulse propagating through the atomic vapor cell exits 
from it earlier than propagating through the same dis- 
tance in vacuum by a time difference that is 315 times of 
the vacuum light propagation time L/c — 0.2 ns. Thus, 
the peak of the pulse exits the cell before it even enters. 
This counterintuitive effect is a consequence of the wave 
nature of light and can be well explained invoking the 
rephasing process in an anomalously dispersive medium. 
The observed superluminal light pulse propagation is not 
at odds with causality, being a direct consequence of clas- 
sical interference between its different frequency compo- 
nents in an anomalous dispersion region. 



through metal is also slower than c, the vacuum speed of 
light. 

These conditions are direct results of the Kramers- 
Kronig relations. Hence, for media in a passive state 
under Eq.(l) where electromagnetic irradiation at all fre- 
quencies are transformed into heat and subsequently dis- 
sipated, anomalous dispersion and transparency cannot 
occur simultaneously in the same frequency region. 

However, for media with gain, the general assumption 
in Eq.(l) no longer holds. In a gain medium such as 
that of an electronic amplifier (e.g. an operational am- 
plifier) or an optical amplifier commonly used for lasers, 
the imaginary part of the electric susceptibility x( v ) can 
become negative over a narrow frequency region where 
gain occurs. Hence, the general results of the inequal- 
ities (2) no longer apply and anomalous dispersion can 
happen in a transparent media, resulting in superluminal 
group velocities that can even become negative. 

A. Classical Theory of Negative Group Velocity 



II. THEORY OF TRANSPARENT ANOMALOUS 
DISPERSION 

For all transparent matter at a thermal equilibrium 
such that it is absorptive at all frequencies in the elec- 
tromagnetic spectrum, the medium's optical dispersion 
is normal |(|. In other words, for transparent media, the 
optical refractive index always increases when the fre- 
quency of the optical excitation increases. Particularly, 
Landau and Lifshitz showed that under the condition 



Im [x( v )] > 0, for any v 



(1) 



and in the special case for media with a magnetic perme- 
ability j.i{v) = 1, two inequalities hold simultaneously: 



. d[n{v)v] 
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> 



n(v) 



(2) 



Here n g is the group velocity index: v g = c/n g , where v g 
is the group velocity. 

In the case of dielectric media, when an incident light 
wave's frequency is below the bandgap where the media 
are transparent, the refractive index n(v) > 1. Hence, 
the group index n g > n > 1, resulting in normal disper- 
sion and the group velocity of light pulses slower than 
c. On the other hand, in the case of metals when the 
incident light's frequency is higher than the plasma fre- 
quency vp, the metal becomes transparent with a refrac- 
tive index n(v) < 1. In this case, the second inequality in 
Eq. (2) becomes more strict and results in n g > 1/n > 1. 
Therefore, the group velocity of a light pulse propagating 



In a transparent medium where absorption is negligi- 
ble, the optical dispersion is primarily governed by its 
refractive index and the dispersion relation can be writ- 
ten: 



k(v) = 2tt 



n{y) ■ v 



(3) 



where k(y) is the wave number for the wave component 
of frequency v. The group velocity is hence given by 



Va = ~ n = Re[ dk ] 



where 



Re[n + vdn/dv] n + vdn/dv' 



dn 
dv 



(4) 



(5) 



is the group velocity index. In all transparent media 
where dispersion is normal, we have dn/dv > 0, result- 
ing in v g < c. In special cases where this normal disper- 
sion is very steep over a narrow frequency region such as 
in the case of Electromagnetically Induced Transparency 
(EIT) ]2q , ^9] that vdn/dv 3> 1, the group velocity can 
be reduced to as slow as 8m • sec -1 |fj0 - 32| L2|| 

Conversely, in a region of transparent anomalous dis- 
persion where vdn/dv <C — 1, a negative group velocity 
index is obtained. If the index of refraction decreases 
rapidly enough with frequency, the group velocity be- 
comes negative. Of course, it is well known that inside an 
absorption line, the refractive index takes a steep drop Q] 
resulting in an "anomalous dispersion" and consequently 
a negative group velocity M. However, the associated 
heavy absorption which results in opacity makes it diffi- 
cult to study these effects. 
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The earlier work of Chiao and coworkers [[7| Jl§|j22[ | pre- 
dicted that it is possible to study steep anomalous disper- 
sion and negative group velocity in a transparent medium 
where gain rather than absorption occurs at the interest- 
ing frequency regions. It is important to have a transpar- 
ent medium where the majority of the electromagnetic 
radiation energy is kept in the light fields rather than 
being dissipated inside the medium via absorption. 

Let us start by considering a classical Lorentz oscillator 
model of the refractive index. The electric displacement 
is given by: 

D = e E + P = eo(l + x)E = E E(l + Na) 1 (6) 

where N is the atomic density and a is the atomic po- 
larizability. The polarization density P = —s^Ne x — 
EoNaE can be obtained using a simple Lorentz model. 

In order to obtain the dipole polarization p — —e x 
for a bound charge with an intrinsic angular frequency 
ujq = Itivq and an angular damping rate T = Attj, we 
start from the equation of motion of the electron: 



x + Tx + uj x = 
Hence, we obtain, 

1 



eE 
m 



eE, 



o & - iu t_ 



eE 



eE 

m oj 2 — loq + iujT 2mujQ ui — u>q + iT/2 ' 



(7) 



(8) 



where the approximation is good as long as wo > T. We 
further obtain for the polarizability, 



1 



1 



a 



2mwo uj — loq + iT/2 Attitiluo 



v — uo + 17 
(9) 



The dielectric susceptibility of the medium thus can be 
written: 



Ne 2 



-/x 



M 



(10) 



where M = Vp/vo with vp being the effective plasma 
frequency and / being the oscillator strength. When two 
absorption lines are placed nearby with equal oscillator 
strengths /i = f 2 — 1, the dielectric susceptibility can be 
written: 



M 



M 



v — V\ 



v — v% 



vy 



(11) 



For a narrow frequency region in the middle between the 
two absorption lines, a steep normal dispersion region oc- 
curs resulting in an ultra-slow group velocity [ ^l| , ^2p^ ] . 

Conversely, for gain lines, a negative oscillator strength 
/ = — 1 is assigned JOl. Hence between two closely 



placed gain lines, the effective dielectric constant can be 
obtained: 



e(v) = \+xiy) = 1 + 



M 



v — v\ + v — v-2 + i"f 



(12) 



For dilute gaseous medium, we obtain from Eq.(12) for 
the refractive index n{v) = n'{v) + in"{v) = 1 + x(i/)/2 
and the real and imaginary parts of the refractive index 
are plotted in Fig.l. It is evident from Fig.l that a steep 
anomalous dispersion region appears without the heavy 
absorption present. In fact, a residual gain persists. Fur- 
thermore, with the correct choice of experimental param- 
eters, the steep drop of refractive index as a function of 
frequency can be made a mostly linear one in this region. 
Thus a light pulse with a frequency bandwidth within 
this narrow linear anomalous dispersion region will expe- 
rience almost no change in pulse shape upon propagating 
through such a medium. 
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FIG. 1. Gain-assisted anomalous dispersion. Figure shows 
frequency-dependent gain coefficient and refractive index. 

Now, we consider the propagation of a light pulse of 
an arbitrary shape through a transparent anomalous dis- 
persing medium of a length L as illustrated in Fig. 2. For 
a scalar light pulse that is decomposed into its positive 
and negative frequency parts: 



E(z,t) = E<- +) (z,t) + E { -\z,t), 
we have for its Fourier decomposition: 



E (+ Hz,t) = 



'2tt 



-i {ujt — k-z) 



(13) 



(14) 



At the entrance of the medium where we denote by z = 0, 
we can rewrite the above expression as 



£ (+) (0,t) 



1 



x e 



2tt 

i (uj— u)o)t 



lUJ0t / d(u- U )E(+\lU - LU Q ) 



(15) 



where u>o is the carrier frequency of the light pulse. Inside 
the transparent anomalous dispersion medium, if over the 
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narrow bandwidth of the incident light pulse E(u> — luq), 
the gain is essentially unity, the propagation is governed 
by the wave vector k(u). Using a Taylor series expansion 
of the wave vector: 



k(oS) — fc(wo) H (w 



. 1 d 2 k 



(uj - UJ ) 2 



(16) 



where the expansion is carried to the quadratic term. 
Nonlinear terms in the expansion of Eq.(16) are often 
associated with "Group Velocity Dispersion (GVD)," or 
chirping terms causing pulse distortion. When these non- 
linear terms in Eq.(16) are negligible, i.e., the dispersion 
is essentially linear, from Eq.(14) and (16) we obtain: 

E {+) (L,t) = _L e -i("ot-fco-L) f d(w-Wo)^ (+) (w-Wo) 
V27T J 

(17) 



x e 



-i (ui—uj )-(t—L/vg) 
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FIG. 2. Pulse propagation through a medium of a length L 
and a group velocity index v g — c/(n + v dn/dv). Pulse prop- 
agation through the same length in vacuum is also shown for 
comparison. 



the vacuum transit time. Hence the output pulse's time- 
varying profile is related to the input pulse by a delay 
that is shorter than the vacuum transit time L/c result- 
ing in a superluminal propagation of the light pulse. 

Furthermore, when the transparent anomalous disper- 
sion becomes stronger to yield: n + v dn/dv = 0, the 
group velocity v g = c/[n + v dn/dv] approaches infinity, 
resulting in a "zero transit time" , such that Eq.(18) gives 
I(L, t) = 1(0, t— L/vg) = 1(0, t). In this case, the output 
pulse and the input pulse vary the same way in time and 
there is no time delay experienced by the pulse propa- 
gating through the medium which has a length of L. 

Finally, when the transparent anomalous dispersion 
becomes very steep such as for the case illustrated in 
Fig.l, the dispersive term v dn/dv which is negative be- 
comes very large in its magnitude such that \vdn/dv\ 3> 
1, resulting in a negative group velocity v g — c/[n + 
v dn/dv] < 0. In this case, Eq.(18) gives I(L,t) = 
1(0, t + \L/v g \), where the quantity \L/v g \ — \n g \L/c 
is positive and can becomes very large compared to the 
vacuum transit time L/c. This means that the intensity 
at the output of the medium of length L, I(L, t), will vary 
in time earlier than that of the input pulse 1(0, t). Thus 
in this case, a "negative transit time" can be observed. 
The time difference between the output pulse and the in- 
put pulse in the form of a pulse advancement, is \n g \ fold 
of the vacuum transit time L/c. Practically, since the 
shape of the pulse is not changed, this results in a rather 
counterintuitive phenomenon where a certain part of the 
light pulse has already exited the medium before the cor- 
responding part of the incident light pulse even enters by 
a time difference that is \n g \ times of the vacuum transit 
time L/c. 

This rather counterintuitive effect is a result of the 
wave nature of light. 



Hence, the intensity of the light pulse as a function of 
time measured with a detector (shown in Fig. 2), I(L, t), 
is related to the incident pulse's time-dependent intensity 
by: 

I(L,t) = £ -^\E^(L,tt = £ -^\E(0,t-L/v g )\ 2 

= I(0,t-L/v g ), (18) 

where A is the beam area. 

Ordinarily, in a normal dispersion medium, the group 
velocity v g < c. Hence, the output intensity of a pulse 
propagating through the medium is retarded by the prop- 
agation time L/v g , resulting in a delay longer than the 
vacuum transit time L/c. In a transparent anomalous 
dispersion medium, the group velocity v g — c/[n + 
v dn/dv] can exceed c provided the anomalous dispersion 
is sufficiently strong such that: n + v dn/dv < 1. In this 
case, the group velocity becomes superluminal: v g > c, 
resulting in a "superluminal transit time:" L/v g < L/c, 



B. Quantum Theory of Atomic Response in 
Transparent Anomalous Dispersion 

In order to correctly model the electromagnetic re- 
sponses of a dilute atomic gaseous medium, one must 
compute its dielectric susceptibility using the quantum 
mechanical treatment. In this section, we provide a sim- 
plified yet realistic model of a transparent anomalous dis- 
persion medium realized using two Raman gain lines in 
a A-system shown in Fig. 3. 

In a simplified A-system, two continuous wave (CW) 
Raman pump light fields are present for the creation of 
the gain doublet shown in Fig.l. The two Raman pump 
fields are different in frequency by a small amount 2Av. 
For simplicity, we will first ignore the Doppler shift and 
assume that the atoms are at rest. 

We begin by treating the simple case of a single Raman 
pump field and compute the linear dielectric susceptibil- 
ity for the Raman probe field. Let us first suppose that 
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all atoms arc initially prepared in an energy ground state 
|1) via optical pumping. In the atomic system, another 
ground state |2) is also present and a Raman transition 
from |1) to 1 2) can take place via an off-resonance two- 
photon transition through an excited state |0) by absorb- 
ing a pump photon and emitting a probe photon, with a 
corresponding transition from the state |1) to |2). The 
effective Hamiltonian for such a A-system coupled with 
the Raman pump and probe fields can be written: 



where 



H = H + H! 



H = -hoj 01 \l)(l\-huj 2\2}(2\, 



(19) 



(20) 



with ujQj = (Eq — Ej)/h for j = 1,2. The interaction 
Hamiltonian can be written: 

Hi = -Tifti e" <Wlt |0)(1| - hQ 2 e-^ 2t |0) (2| + H.c, (21) 

where VLj = (e 3 • fj,oj)Ej/2U (j — 1,2) is the Rabi fre- 
quency for the pump and the probe fields, respectively. 
jiQj is the dipole moment between the states \j) and the 
excited state |0). Because Ex is a strong Raman pump 
field while E 2 is a weak Raman probe field, therefore, we 
have > |fi 2 |- 



F = 4,m = -3) 



6P 3/2 |0) 




F = 4,m = -2) 



|F = 4,m = -4) 



FIG. 3. Schematic atomic level diagram. 

Without involving complicated density operator treat- 
ment which is only necessary in a full quantum theory 
where the electromagnetic field is also quantized, we can 
apply here a simplified state vector treatment in order 
to obtain the dielectric susceptibility of a dilute gaseous 
atomic medium. The state vector of an atom can be 
written: 

\m) = ao(t) |0) + ai(i) e^ 0lt |1) + a 2 (t) e^ 2 * |2). (22) 

Hence, we obtain for the amplitude of the excited state 
|0>, 



a (t) = i fiie _i Alt ai + i Q 2 e' lA2t a 2 , 



(23) 



where Aj = uoj — u>oj, (j = 1, 2) represent the detuning 
of the Raman pump and probe beams, respectively. To 
the lowest order of approximation, we have for the state 
amplitudes a\ w 1 and a 2 ~ 0. We further note the fact 
that the common detuning Ao = (Ai + A 2 )/2 is much 
greater than the differential detuning Ai — A 2 . Under 
these conditions, we have from Eq.(23): 



a (t) 



■^e- iA ^ai. 
Aq 



(24) 



Note here that we have omitted the decay of the excited 
state amplitude ao in Eq.(23). First, the broadened line 
width is still much narrower than the common detuning 
A such that Eq.(24) is a good approximation. Second, 
in a Raman scheme as used in the experiment, the ma- 
jor factor of the Raman transition broadening is due to 
transit time broadening rather than the effective excited 
state decay rate: (ilf/ Aq)/T 2 , where T 2 is the excited 
state lifetime. Hence, Eq.(24) is a good approximation. 

Next, we obtain the equation of motion for the ampli- 
tude of the Raman final state |2): 



62(f) 



-t-^2 e -*( A i-*')*Oi 



Ao 



r a 2 . 



(25) 



Here, a phenomenological decay rate T is used to ac- 
count for the Raman transition line broadening. Solving 
Eq.(25), we obtain 



02 (t) = 



1 



Ao J (Ai - A 2 



-i(Ai-A 2 )t 



iT 



a x . (26) 



Using Eq.(24) and Eq.(26) and the definition of the di- 
electric polarization P = N /i 2 ,o P02 = X £ o E 2 where the 
density matrix element P02 = ao a 2 e ~ J " 02t , we obtain for 
the dielectric susceptibility for the Raman probe field E 2 , 



X(A 2 ) = -JV- 



IM02 ■ e 2 | 
2he a 



A 2 
^0 



1 



(Ai- A 2 )-*r 



(27) 



Using the fact that Ai — A 2 = — 2tt \v 2 — (yx —^01 + ^02)], 
and redefining v = {vx — vqx + ^02) while replacing v 2 
with u, we can rewrite Eq.(27) as: 



M 



v - i/ +«7 



where the factor 



M = N ■ 



1^02 • e 2 | 2 |Oi| 
A-kTlex) Al 



(28) 



(29) 



When two Raman pump beams of equal strengths and 
slightly different frequencies v§ — Av and vq + Av are 
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present, the dielectric susceptibility of the Raman probe 
field becomes: 



M 



M 



{v — i>q — Ai/) + i 7 {y — vq + Aj/) + 17 



, (30) 



similar to that obtained using a classical Lorentz model 
given in Eq.(ll). 

However, in a gaseous atomic medium, there is Doppler 
broadening that demands the expression in Eq.(29) and 
Eq.(30) being modified. Specifically, the common detun- 
ing in Eq.(29), Ao, needs to be replaced by the shifted 
detuning Ao + vq V/c, where vq is the carrier frequency 
of the light field and V is the velocity of those atoms 
in a specific velocity group G(V) for which the atoms 
move along the light propagation direction with a speed 
V. Owing to the co-linear propagation geometry, the 
Doppler shifts in both Ai and A2 are essentially the 
same. Hence, the shift in frequencies in Ai — A2 or 
v — vq zb Ais are essential zeo to within an error of Aj/ • V/c 
which is negligible. Therefore, Eq.(30) need only be mod- 
ified by replacing the expression of the M-factor with: 



M = N ■ 



IM02 • e 2 [ 



dV 



(A + v ■ V/c)- 



■G{V). (31) 



The effects of this modification are two fold. First, the 
quadratic dependence on 1/(A + v ■ V/c) 2 in Eq.(31) is 
an even function which prevents cancellation of the ef- 
fect due to Doppler broadening. This is a direct result 
of using Raman transitions. Second, for some velocity 
groups where the shifted detuning Ao + v ■ V/c may van- 
ish. However, for these velocity groups, the Raman pump 
beams act like reversed optical pumping which empty 
these velocity groups such that G(Vo) — 0. This is similar 
to the "spectral hole burning" effect commonly know in 
laser physics. Furthermore, the atoms reversely pumped 
away in these velocity groups act like a broadband weak 
absorber that helps to compensate the residual gain as 
shown in Fig.l. 

The spectral linewidths of the gain lines are approxi- 
mately given by 
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vr 
2ir wo ' 



(32) 



due to transit broadening |33| ]. Here vr is the mean 
atomic velocity in the radial direction and wq is the mean 
radius of the Gaussian beam. 



C. Parameter Dependence in Gaussian Pulse 
Propagation 

The primary result of Eq.(18), I(L,t) = I(0,t-L/v g ), 
applies to light pulses of arbitrary shapes, provided that 
over the pulse bandwidths the dispersion is transparent 
and linear. In order to further examine the implications 



and limits of pulse propagation in a transparent anoma- 
lous dispersion medium, here we treat the propagation of 
a Gaussian pulse. 

For the propagation of a Gaussian pulse through a 
transparent anomalous dispersion medium, there are 
three important parameters: the peak intensity, the pulse 
advancement, and the pulse shape distortion. In this 
section, we calculate the parameter dependence of these 
factors on various experimental parameters. These pa- 
rameter dependence will be compared with experimental 
results. 

To begin, let us refer to Fig. 2 and consider a Gaussian 
pulse of a temporal duration r and an angular carrier 
frequency uiq, for which the incident waveform becomes: 



E(z, t) = E e 



-{t-z/c) 2 /2T 2 e -iu> (t-z/c) 



[z < 0). (33) 



At the entrance of a transparent anomalous dispersion 
medium z — 0, the electric field hence varies in time as: 



E(0,t) = E Q e- t2/2r2 



-ILJQt 



(34) 



Therefore, at the exit surface of a medium of length L, 
the electric field varies in time: 



where 



tE e- ,lL 



-T 2 {u)-u> ) 2 /2 e -i(u t-k I^fa 



a Eqc 



-0 2 (t-L/v g ) 2 /2T 2 -!u„ [t-n'{wo)L/c] 



(35) 



a = e -n"(u )u L/c-r,L . ^ = &g L . ^ 

describes the change in pulse amplitude due to residual 
amplification and a weak broadband absorption factor 77, 
resulting in a phenomenological gain coefficient gL. The 
change in the width of the pulse is entirely characterized 
by the parameter 



L d 2 k 

T 2 dhJ 2 



(37) 



Using the simple expression in Eq.(30) and (31), here 
we calculate the scaling dependence of various quantities 
which can be compared with experimental results. 

First, using Eq.(30), we obtain a residual gain coeffi- 
cient of 



-n"(u> ) luq/c- rj 



2ir vq M 



7 



(Azy 2 + 7 2 ) 



17- 



(38) 



We note here that the factor M is linear in the atomic 
density N and the Raman pump beam power |f2i| 2 as 
given in Eq.(31). Meanwhile, the residual loss coefficient 
?7, which is due to absorption of atoms reversely pumped 
away into other atomic states is also linear in the atomic 
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density N. Its dependence on the Raman pump power 
is more complex since power broadening and common 
detuning are involved. 

Next, we calculate the pulse advance (negative delay), 
AT, for a Gaussian pulse propagating through an anoma- 
lous dispersion medium of length L. Referring the ex- 
perimental situation illustrated in Fig. 2, a light pulse 
propagating through a transparent anomalous dispersion 
medium of a length L is advanced compared to the same 
pulse propagating through the same distance in a vacuum 
by, 



L dn 

AT = (1 - rig) — = -p ~r 

c dv 



L 

c 



(39) 



Using the expression in Eq.(30), we have for the pulse 
advance: 



L Av 2 - 7 2 

AT = U M — ■ — ^ 7T-T . 

c (A^ 2 +7 2 ) 2 



(40) 



The pulse propagating through the transparent anoma- 
lous dispersion medium is advanced provided that Av > 
7. Here we further note that the pulse advance, AT, is 
linear in M which is linear in the atomic density N and 
the Raman pump power |r2i| 2 . These parameter depen- 
dences will be compared with experimental results. 

Finally, we compute the pulse width change factor /?. 
We notice from Eq.(37) that j3 depends on the derivative 
of the group velocity index n g . From Fig.l, it is appar- 
ent that the frequency-dependent change of the real part 
of the refractive index n (v) is essentially linear in the 
narrow frequency region between the two gain lines. It is 
the imaginary part, i.e., the gain coefficient that appears 
to have a quadratic change as a function of frequency. 
Therefore, the majority of the potential change in the 
width of an incident pulse will be due to the quadratic 
change in the imaginary part of the refractive index. We 
obtain using Eq.(30), 



d 2 k 
du 2 



= -iMv 



7 (3A^ 2 -7 2 ) 
irc(Av 2 + 7 2 ) 3 



(41) 



Hence, from Eq.(7), the pulse width distortion factor f3 
becomes: 







1 - Mv 



L 7 (3Av 2 - 7 2 



7TC T 2 (A^ 2 +7 2 ) 3 



(42) 



When various parameters such as the M factor (linear in 
Raman pump power and atomic density) and the length 
of the atomic chamber L, are significantly increased, or 
the pulse duration r is significantly decreased, the fac- 
tor [3 will become reasonably different from (larger than) 
unity. In these cases, f3 > 1 will result in a narrowing in 
the pulse's temporal width. 

In the present experiment, this pulse distortion factor 
(3 is only 1.002. 



III. EXPERIMENTS 



From the theoretical treatment above, we see that a 
variety of parameters have to be taken into considera- 
tion in experimentation. From an experimental point of 
view, one must satisfy the following requirements. First, 
a gain doublet must be obtained for which the anoma- 
lous dispersion between the gain lines can become linear 
to avoid any pulse distortion. In previous work, excited 
state population inversion was considered to obtain gain 
p2| . However, spontaneous emission and the short ex- 
cited state life time would cause such gain doublets to be 
very difficult to sustain. It is important to have a steady 
state gain with a lifetime longer than the pulse duration 
t to avoid transient effects and the associated compli- 
cations. Second, the medium must be transparent since 
opaque anomalous dispersion has been long known and 
has resulted in controversies in terms of interpretations. 
While ideally the dispersion shown in Fig.l is transpar- 
ent, residual absorption and the associated loss are often 
present and cannot be simply discarded in experimental 
situations. Third, in order to show superluminal light 
pulse propagation in a linear regime, one must employ 
a very weak light pulse for which the photon number is 
far less than the atomic number in order to avoid Raman 
gain saturation. Fourth, in order to achieve a reason- 
able accuracy in the measurement, a system should be 
designed to demonstrate a negative group velocity. In 
this case, the pulse advancement under conventional ex- 
perimental situation will be substantially large compared 
with commonly obtained accuracy (about Ins). A num- 
ber of other experimental conditions also have to be con- 
sidered such as atomic density, polarization decay time, 
etc. and they will be discussed in the following sections 
as well. 



A. Experimental Setup 



The experiment is performed using an atomic Cesium 
(Cs) vapor cell at 30°C and the main setup is shown in 
Fig. 4. The cesium atoms are confined in a 6-cm-long 
Pyrex glass cell coated with paraffin for the purpose of 
maintaining atomic ground-state spin polarization. The 
atomic cell is placed inside a magnetic shield made of 
a thin layer of high-// metal material inside which the 
Earth magnetic field is reduced to sub-milli-Gauss level. 
A Helmholtz coil (not shown in Fig. 4) produces a uni- 
form magnetic field inside the magnetic shield parallel 
to the propagation direction of all optical beams. This 
uniform field is approximately 1 Gauss serving the pur- 
pose of defining a quantization axis for optical pumping. 
Inside the magnetic shield, the air temperature is con- 
trolled using a heater servo system in order to control 
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the temperature of the Cesium cell. During data acquisi- 
tion, this control system is turned off to avoid any stray 
magnetic field. Having good thermal insulation, the tem- 
perature of the atomic cell remains the same during the 
data acquisition time. 

In region-I of Fig. 4, two optical pumping laser beams 
prepare almost all Cesium atoms into the ground-state 
hyperfine magnetic sublevel 6S1/2, \F — 4, m = —4) that 
serves as the state |1) in Fig. 3. Laser-1 is a narrow 
linewidth diode laser locked to the 852-nm D2 transi- 
tion of Cs using a Lamb-dip technique and empties the 
6S1/2, F = 3 hyperfine ground statesLaser-2 is a broad- 
band tunable Tksapphire laser tuned to the 894-nm D\ 
transition of Cesium. The linewidth of laser-2 covers 
transitions from both the 651/2, F — 4 and F = 3 hyper- 
fine ground states to the 6P1/2 excited state. Both laser 
beams are initially linear polarized and are turned into 
left-hand polarization (a— ) using a quarter- wave plate 
placed before the atomic cell. Inside the vapor cell, Ce- 
sium atoms collide with the paraffin-coated glass walls, 
the atoms change their velocities inside the Doppler pro- 
file. However, their ground state spin polarizations are 
not changed during collisions. Hence, all atoms inside the 
entire Doppler broadening profile are optically pumped 
into the ground state \F = 4, m = —4) quickly. The 
mean dephasing time of the ground state spin polariza- 
tion of Cesium atoms in a paraffin coated cell is of order 
a fraction of a second. 
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FIG. 4. Schematic experimental set up. Two optical pump- 
ing beams tuned to the Cesium (Cs) atomic Dl and D2 transi- 
tions prepare the atoms in its ground state hyperfine sublevel 
\F — 4, m — —4). Two Raman pump beams and a Raman 
probe beam derived from a common narrow linewidth diode 
laser propagate colinearly parallel to a small magnetic field B 
through the atomic cell. Two A/4-plates (QWP1 and 2) are 
used to prepare the three light beams into the correspond- 
ing circular polarization states and then separate them for 
analysis. 

In region-II, three light beams derived from the 
same narrow linewidth diode laser propagate colinearly 
through the cell. Two strong continuous-wave (CW) Ra- 



man pump beams are right-hand polarized (er+) and a 
weak Raman probe beam is left-hand polarized (a— ). Us- 
ing 3 acousto-optical modulators (AOM's), the frequency 
difference of the two Raman pump beams can be tuned 
continuously over a few MHz while the probe beam can 
also be tuned in frequency and can be operated in both 
CW or pulsed mode. The typical carrier frequency of 
the AOM's is 80 MHz and the linewidth is about 20 kHz. 
A residual optical beam that is shifted in frequency by 
80 MHz generated from the same AOM which modulates 
the probe beam is also available for the refractive index 
measurement. 



B. Experimental Methods and Results 

First, we operate the Raman probe beam in a tunable 
CW mode to measure the gain and refractive index of 
the atomic system as a function of the probe frequency 
detuning. Fig. 5 shows the measured gain coefficient and 
the refractive index. In order to obtain the gain coeffi- 
cient, we first measure the intensity of the transmitted 
probe beam as a function of probe frequency. We then 
extract the gain coefficient. 
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FIG. 5. Measured refractive index and gain coefficient. 
The superposed curve over the index data is obtained using 
Eq. (30) with parameters i>o, Av, and 7 obtained experimen- 
tally. The length of the cell is 6 cm. 

The refractive index is measured using a radio- 
frequency (rf) interferometric technique. We send two 
CW light beams from the probe beam input port of the 
polarized beamsplitter. These two optical beams are dif- 
ferent in frequency by Q w 80 MHz (derived from the 
zeroth and the first order output of the probe AOM). A 
weak probe beam is tuned in frequency to be close to 
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the Raman pump frequency and a strong local oscillator 
field shifted by £1 s=s 80 MHz is outside of the interesting 
frequency range. The rf beating signal is detected us- 
ing two fast detectors Dl and D2. We AC couple both 
output signals and record the beating signal of D2 

V on {t,cj) on (v)) = V cos[flSt - cj) on {v) - 4>q] (43) 

as a function of the AOM frequency f2. using the AC 
coupled Dl voltage signal as the trigger. We then apply 
a least-squares fitting procedure using Eq.(43) to obtain 
the frequency dependent phase shift 4> on (v). In Eq.(43), 
St is a residual difference between the electrical delays of 
the trigger channel and the signal channel and is mini- 
mized. cf>o is a fixed phase factor. Of course, the fitting 
procedure always only yields the combined phase shift, 



t (u) -6(1 St = 2im{v) L/X - SQ St 



(44) 



where S(t = (I — (Iq (^o=80 MHz) is the change in fre- 
quency relative to the mean carrier frequency. In order 
to achieve high accuracy of the refractive index measure- 
ment, one must separately measure the factor 8(l8t. To 
do this, wc tune the diode laser outside of the Doppler 
profile where the atoms are irrelevant and Eq.(43) be- 
comes: 



V off(t) = Vi cos[(lSt - (p off - cj> ] 



(45) 



since the phase shift here <j) ff = 0. Hence, we di- 
rectly measure the factor SClSt and can subtract it out 
from the combined phase shift in Eq.(44) to obtain the 
the refractive phase shift to obtain the refractive index: 
<t>oniy) = 2-kti(v) L/X. The result is shown in Fig. 5. The 
superimposed curve is obtained from Eq.(30) using pa- 
rameters obtained from the gain measurement. From 
Fig. 5, we readily see that a negative change of An= - 
1.8xl0~ 6 in the refractive index occurs over a narrow 
probe frequency range of IS.v~ 1.9 MHz between the two 
gain lines. Using the expression of the group-velocity 
index, we obtain the result n g = —330 (±30) in that fre- 
quency region. The 10% error reflects the accuracy of 
the phase measurement. 

Next, a pulsed Raman probe beam is employed to ob- 
serve the superluminal propagation. A near Gaussian 
probe pulse with a 2.4 /j,sec FWHM is generated by ap- 
plying a biased sinusoidal electronic signal to the probe 
beam A/O modulator. The repetition rate is 50 kHz. A 
portion of the pulsed probe beam is divided at a beam- 
splitter before the atomic cell and aligned onto photo- 
diode Dl as a reference. Because the total number of 
atoms in the probe volume limits the maximum net en- 
ergy gain of the probe pulse, we use a very weak probe 
beam (w 1//W) in order to avoid saturation and hence 
to optimize the anomalous dispersion. A high sensitivity 
avalanche photodiode, reverse-biased below breakdown, 
serves as detector D2 to measure the weak probe pulse 



that propagates through the atomic cell. The photoelec- 
tric current produced by detector D2 is converted to a 
voltage signal using a 500-J1 load resistor and recorded 
by a digitizing oscilloscope using a synchronized output 
signal from the pulse generator as the trigger. Pulses 
from detector Dl are also recorded. 

In order to measure the pulse propagation time, we 
first tune the diode laser that produces the Raman pump 
and probe beams far off-resonance from the 852 nm 
Cesium D2 lines (by 2.5 GHz) to measure the time- 
dependent probe-pulse intensity. When the laser is 
placed far off-resonance, the atoms have no effect and 
the probe pulse propagates at the speed c inside the cell. 
We then tune the diode laser back to within the Doppler 
absorption profile and lock it on its side. Using the same 
synchronized pulse generator output signal as the trig- 
ger, we record the time-dependent probe pulse intensity 
measured by detector D2. We verify that no systematic 
drift is present by tuning the laser off-resonance again by 
the same amount and record the probe pulse signal; the 
two off-resonance pulses are identical to within less than 
1 ns. Probe pulses both on and off-resonance are shown 
in Fig. 6 (average of approximately 1000 pulses). 



-1.4 -1.2 -1.0 




FIG. 6. Measured pulse advancement for a light pulse 
traversing through the Cesium vapor. (A), light pulse far 
off-resonance from the Cesium D2 transitions propagating at 
speed c through 6 cm of vacuum. (B), same light pulse prop- 
agating through the same Cs-cell near resonance with a neg- 
ative group velocity — c/315. Insets show the front and trail 
parts of the pulses. A and B are both the average of 1000 
pulses. Off-resonance pulse (A) is normalized to the magni- 
tude of (B). 

Curve A in Fig. 6 shows the light pulses far off- 
resonance from the Cesium D2 transitions propagating 
at speed c through 6 cm of vacuum. Curve B shows 
the same light pulse propagating through the same Cs- 
cell near resonance at a negative group velocity —c/315. 
Probe pulses on resonance show a 40% transmittancc and 
this is due to the broadband absorption of those atoms re- 
verse pumped away from the state \F = 4, m = —4). It is 
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evident that there is almost no change in the pulse shape. 
The front edges and the trailing edges of the pulses are 
shown in the insets; both edges are shifted forward by 
the same amount. 

Using a least square fitting procedure, we obtain a 
pulse advancement shift of 63 (±1) nsec. Compared 
with the 0.2 nsec propagation time for light to tra- 
verse the 6-cm length of the atomic cell in vacuum, the 
63 nsec advancement gives an effective group index of 
rig = 315(±5). This is in good agreement with that in- 
ferred from the refractive index measurement. The pulses 
measured with detector Dl are also recorded in the se- 
quence of the off-, on-, off-resonance pulse propagation 
measurements and are found to be identical to within 
1.5 ns. 

We note here that the measured supcrluminal pulse 
propagation inside the transparent anomalous dispersion 
medium is a linear effect. We further estimate the pho- 
ton number per pulse and the interacting atomic number 
to show that there is no saturation effect present. In the 
experiment, the measured voltage signal peak strength is 
V p = a£ RG Nphhtoo, where a w 0.5A/W is the photo 
responsivity of the avalanche photo detector and £ = 0.2 
is an effective efficiency of the detection imaging system. 
R=500 is the load resistance. G w 80 is the avalanche 
gain, hcdo = l.beV is the photon energy. Hence, we ob- 
tain the peak photon rate N p h «5x 10 12 / sec. Using the 
2.4 [isec FWHM as the pulse duration, each probe pulse 
contains approximately 1.2 x 10 7 photons. On the other 
hand, inside the volume of the probe light beam ttwq L, 
there are on average N ttwq L atoms at any given mo- 
ment, where N s=s 10 n cra -3 is the atomic density. The 
beam radius is approximately 90 /im and the Cesium cell 
is of a length 6 cm, the atoms inside the beam volume is 
approximately 1.4 x 10 8 at any given moment. However, 
since atoms are coming in and out of this volume within 
an average time of 2 wq/Vr rs 1 /isec, during the 2.4 fisec 
pulse duration, there are approximately 3.4 x 10 8 atoms 
inside the volume of the light pulse, much larger than 
the photon number per pulse. Hence, gain saturation 
effects are insignificant and the observed superluminal 
pulse propagation is a linear effect. 

To further analyze the linearity of the response of the 
anomalous dispersion medium, we compute the pulse 
area of the Raman transition ■ r. Here = 
IO1O2/A0I is the effective Raman Rabi frequency. By 
using Eq.(29) and by noting that ^2 = (M02 ■ e 2 )E 2 /2h, 
we obtain 



Hence we obtain from Eq.(46) 



M\E 2 \ 



N ■ 



TT£ 



A 



(46) 



where N is the atomic density. Furthermore, we have for 
the probe electric field E 2 , P 2 — 2e a cA\E 2 \ 2 — N p h Tuoq, 
where P 2 is the Raman probe power, A = nw 2 is the 
beam cross-section, and Tiloq is the energy per photon. 



irMuj N ph 
2NcA 



We further recall from Eq. (28) that the peak Raman gain 
G = gL at one of the resonances is G = ttML/jX. 
Therefore, we obtain for the effective Raman Rabi fre- 
quency 



9.1 



Ao 



/ TrG-/N ph 
NitWqL ' 



(47) 



Using the experimental parameters of N p h — 5 x 
10 12 /sec, NitwqL = 1.4 x 10 8 , peak gain G = 0.7, 
and 7 = 0.45 MHz, we obtain Vl R = 1.25 x 10 5 /sec. 
This is far smaller than the typical Raman detuning 
IS.v w lAMHz. Recall that during the 2.4 /zsec pulse 
duration, the average dwell time for the atoms is At = 
2wo/Vr = 1/xsec. Hence, the equivalent pulse area be- 
comes riflAr Rj 0.13. It is apparent that saturation of 
Raman gain is negligible. 

Next, we perform experiments to test the parameter 
dependence given by Eq.(40). Namely, the pulse ad- 
vancement I AT I is linear in the parameter M which de- 
pends linearly on the Raman pump power Pr and the 
atomic density TV. First, the pulse advancement AT is 
measured as a function of the total Raman pump power 
of the two gain lines using the method described above. 
Fig. 7 shows the experimental result. In Fig. 8, we show 
the measured parameter dependence of pulse advance- 
ment I AT I as a function of atomic density. We measure 
the optical density of the atomic medium in an ancillary 
linear absorption experiment while varying the tempera- 
ture of the atomic cell. The optical density of the atomic 
cell is proportional to the atomic density N. 



100 




200 400 600 

Total Raman Pump Power (jiW) 



800 



FIG. 7. Measured pulse advance |AT| dependence on total 
Raman pump power. Solid line shows the least-square fitting 
based on the linear dependence shown in Eq. (40) 
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FIG. 8. Measured pulse advance |AT| dependence on the 
optical density of the Cs chamber (linear in atomic density 
TV). Solid line shows the least-square fitting based on the 
linear dependence shown in Eq.(40) 
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FIG. 9. Measured pulse gain coefficient and pulse de- 
lay/advance AT as a function of the carrier frequency of the 
light pulse. Dashed line shows the theoretical curve of refrac- 
tive index using parameters extracted from the gain coefficient 
using Eq.(30). Solid line shows the theoretical value of pulse 
advance/delay based on the gain coefficient and refractive in- 
dex. 

Finally, we measure the pulse delay/ advancement AT 
as a function of the pulse carrier frequency v. Fig. 9 shows 
the measurement results. For several values of the car- 
rier frequency the pulse gain is measured. We further 
record the pulse forms for the laser frequency tuned off- 
and on- resonance. Using the measured pulse gain, we 
obtain the gain coefficient and fit the data using Eq.(30). 
Using the fitting parameters, we compute the theoretical 
value of refractive index and group delay/advancement 
plotted also in Fig. 9. The large error bar for the pulse 
delay measurement near the gain lines are due to the 
rapid change in group velocity index and gain, causing 
severe pulse shape distortions. The simple model based 



on a Lorentz oscillator is in good agreement with the 
experimental results. 



IV. DISCUSSION AND CONCLUSION 

As remarked by Lord Raylcigh [Q, the group veloc- 
ity of a light pulse is the result of interference between 
its various frequency components. Here we note that 
the measured negative and superluminal group veloc- 
ity of a light pulse propagating through a transparent 
anomalous dispersion medium is due to the physical effect 
of "rephasing." Specifically, inside an anomalous disper- 
sion medium, a longer wavelength (redder) component 
of a light pulse has a slower phase velocity, contrary to 
the case of a normal dispersion medium. Conversely, a 
shorter wavelength component (bluer) has a faster phase 
velocity. Inside a medium of refractive index n, the ef- 
fective wavelength of a light ray is modified: A' = A/n, 
where A is the vacuum wavelength. Therefore, in a suffi- 
ciently strong anomalous dispersion medium, the redder 
incident ray will have a shorter wavelength and hence be- 
comes a bluer ray, while an incident bluer ray will have 
a longer wavelength to become a redder ray. This results 
in an unusual situation where the phases of the different 
frequency components of a pulse become aligned at the 
exit surface of the medium earlier than even in the case 
of the same pulse propagating through the same distance 
in a vacuum. 

This highly unusual circumstance was poorly under- 
stood because it only occured in absorption lines previ- 
ously. Furthermore, inside an opaque medium, the ma- 
jority of the incident light energy is absorbed and sub- 
sequently dissipated by the medium making it difficult 
to define the energy velocity of light. In the present ex- 
periment, an anomalous dispersion region is created in a 
transparent medium. Hence it becomes possible to speak 
about the energy velocity of a light pulse J|4|,|35|] . 

Here we note that the physical mechanism that governs 
the observed superluminal light propagation has been 
traditionally viewed as "virtual" reshaping |t],|25| . In the 
past, such superluminal light pulse propagation has been 
widely viewed as the result of the amplification of the 
pulse's front edge and the absorption of its tail, inspite it 
had been repeatedly pointed out that such reshaping is 
actually a "virtual process" Jt],^5| . In the present exper- 
iment, the 2.4/isec (FWHM) probe pulse has only a 160 
kHz bandwidth (FWHM) that is much narrower than the 
2.7 MHz separation of the two gain lines and the probe 
pulse is placed in the middle of these gain lines spectrally. 
Hence, the probe pulse contains essentially no spectral 
components that are resonant with the Raman gain lines 
to be amplified. Therefore, the argument that the probe 
pulse is advanced by amplification of its front edge does 
not apply. Furthermore, the average time an atom stays 
inside the volume of the Raman probe beam is < 1/xsec, 
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shorter than the 2.4/isec FWHM of the pulse. Hence, 
even if the atoms are amplifying the probe pulse, both 
the leading and the trailing edges would be amplified as 
both edges interact with atoms in the same state. This 
is not consistent with the experiment. Hence it is worth 
emphasizing that the "reshaping" of the pulse can be only 
viewed as "virtual reshaping" |7],|25| aimed at providing 
an intuitive understanding. Strictly speaking, the super- 
luminal light propagation observed here is the result only 
of the anomalous dispersion region created with the as- 
sistance of two nearby Raman gain resonances. In other 
words, it is largely due to the rapid anomalous change of 
the refractive index as is shown in Eq.(39) rather than 
the gain. When the gain coefficient becomes large, its ef- 
fect appears as the compression of the pulse, as indicated 
by Eq.(37) and (40). We further stress that the observed 
superluminal light propagation is a result of the wave 
nature of light || . It can be understood using the clas- 
sical theory of wave propagation in an anomalous dis- 
persion region where interference between different fre- 
quency components produces this rather counterintuitive 
effect. 
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